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Abstract 


A method  is  developed  for  defining  localized  states  and  effective 
Hamiltonians  in  perturbed  crystals.  The  method  is  based  on  the 
localization  ideas  in  the  kq  representation  for  perfect  lattices. 

An  equation  is  derived  defining  localized  states  for  perturbations 
caused  by  an  impurity,  the  magnetic  and  electric  fields.  First, 
the  impurity  problem  is  considered  in  detail.  A correction  term 
is  obtained  to  the  one-band  Koster-Slater  effective  Hamiltonian  . 

It  is  shown  to  be  significant  for  bound  states  and  scattering  cross  1 
sections  of  a localized  impurity.  Second,  an  orthonoraal  set  of 
localized  states  for  a crystal  with  a perturbation  is  developed. 

It  includes  the  impurity  problem,  surface  states,  super lattices  and 
other  perturbations.  These  localized  states  are  used  for  deriving 
one-band  effective  Hamiltonians  up  to  second  order  of  the  perturbation 
expansion.  Relatively  simple  results  for  localized  states  and  one- 
band  Hamiltonians  are  obtained  in  the  cases  of  wide  and  narrow  energy 
gap  crystals.  The  orthonormal  set  of  states  that  diagonalize  the 
Hamiltonian  is  also  used  for  deriving  an  expression  for  local  charge 
densities  in  a perturbed  cyrstal  which  can  directly  be  compared  with 
experiment. 

Third,  the  symmetric  coordinates  in  solids  are  used  for  developing  Wannier 
functions  in  the  presence  of  a magnetic  field.  The  functions  are  shown 
to  form  an  orthonormal  set  with  respect  to  different  bands  and  different 
sites  of  the  crystal.  A simple  relation  is  established  between  the 
eigenfunctions  of  a Bloch  electron  in  a magnetic  field  and  the  newly 
developed  Wannier  functions.  This  relation  is  used  for  an  entirely 
elementary  derivation  of  one-band  effective  Hamiltonians.  Multiband 
operators  are  defined  in  the  space  of  envelope  functions  for  different 
physical  quantities,  e.g.  the  radius  vector,  velocity  and  acceleration. 
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1.  ABSTRACT 

A method  is  developed  for  defining  localised  states  and 
effective  Hamiltonians  in  perturbed  crystals.  The  method  is  based 
on  the  localization  ideas  in  the  kq  representation  for  perfect 
lattices.  An  equation  is  derived  defining  localized  states  for  per- 
turbations caused  by  am  impurity,  a magnetic  and  electric  fields. 
First,  the  impurity  problem  is  considered  in  detail.  A correction 
term  is  obtained  to  the  one-band  Koster-Slater  effective  Hamiltonian. 
It  is  shown  to  be  significant  for  bound  states  and  scattering  cross 
sections  of  a localized  impurity.  Second,  an  orthonormal  set  of 
localized  states  for  a crystal  with  a perturbation  is  developed.  It 
includes  the  Impurity  problem,  surface  states,  superlattices  and 
other  perturbations.  These  localized  states  are  used  for  deriving 
one-band  effective  Hamiltonians  up  to  second  order  of  the  perturba- 
tion expansion.  Relatively  simple  results  for  localized  states  and 
one-band  Hamiltonians  are  obtained  in  the  cases  of  wide  and  narrow 
energy  gap  crystals.  The  orthonormal  set  of  states  that  diagonalize 
the  Hamiltonian  is  also  used  for  deriving  an  expression  for  local 
charge  densities  in  a perturbed  crystal  which  can  directly  be  compared 
with  experiment. 

Third,  the  symmetric  coordinates  in  solids  are  used  for 
developing  Wannler  functions  in  the  presence  of  a magnetic  field. 

The  functions  are  shown  to  form  an  orthonormal  set  with  respect  to 
different  bands  and  different  sites  of  the  crystal.  A simple  relation 
is  established  between  the  eigenfunctions  of  a Bloch  electron  in  a 
magnetic  field  and  the  newly  developed  Wannier  functions.  This 
relation  is  used  for  an  entirely  elementary  derivation  of  one-band 
effective  Hamiltonians.  Milt ib and  operators  are  defined  in  the 
space  of  envelope  functions  for  different  physical  quantities,  e.g. 
the  rsdius  vector,  velocity  and  acceleration. 


i 


2. 


INTRODUCTION 


The  recent  interest  in  the  Koster-Slater  impurity  problem^" 

2 

was  caused  by  a new  idea  put  forward  in  a paper  by  Kohn  and  Onffroy  . 
The  idea  consists  in  defining  one-band  generalized  Wannier  functions 
(CWF)  that  are  suitable  for  expanding  the  corresponding  one-band 
eigenfunctions  of  a perturbed  crystal.  By  doing  so  the  multiband 
Koster-Slater  equation  reduces  to  a one-band  problem.  In  addition, 
the  local  density  of  states  and  the  charge  density  can  directly  be 
expressed  in  the  GWF. 

In  this  work  an  equation  for  defining  localized  states  and 

effective  Hamiltonians  in  perturbed  crystals  is  developed.  The 

equation  is  written  in  the  framework  of  localized  states  in  perfect 

crystals  in  the  kq  representation.  The  localization  problem  on 
3 4 

ideal  lattices  * was  recently  solved  and  the  approach  here  will  be 
an  extension  of  Ref.  4 to  perturbed  crystals.  A variety  of  perturba- 
tions is  considered,  e.g.  the  impurity  problem,  the  magnetic  and 
electric  fields.  In  developing  the  above-mentioned  equation  some 
simple  facts  about  the  localization  problem  in  the  kq-representation 
were  used.  First,  a localized  function  in  a perfect  crystal,  a Wannier 
function,  for  some  band  in  the  kq-representation  is  equal  to  the 
Bloch  function  for  the  same  band.  Secondly,  localized  functions  in  the 
kq-representation  satisfy  the  Bloch  equation  for  a perfect  crystal. 
Thirdly,  Wannier  functions  in  the  kq-representation  for  different 
sites  differ  from  one  another  by  a simple  exponential  factor.  In  a 
non-perfect  crystal  these  facts  can  easily  be  generalized,  and  one 
is  then  lead  to  an  equation  for  perturbed  localized  states.  Because 
of  their  localized  character,  one  expects  the  perturbed  localized 
states  not  to  differ  very  strongly  from  the  corresponding  localized 
states  in  a perfect  crystal.  Having  this  in  mind,  a pertubatlve 
procedure  was  developed  in  for  deriving  localized  states  and 
effective  Hamiltonians  in  perturbed  crystals. 
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Using  the  set  of  localized  orthonormal  functions  a general 
theory  is  developed  for  deriving  one-band  effective  Hamiltonians. 

The  results  are  obtained  up  to  second  order  in  the  perturbation 
expansion.  The  orthonormal  functions  that  lead  to  the  one-band 
effective  Hamiltonians  are  shown  to  be  suitable  for  deriving  local 
charge  densities  in  a crystal  with  a perturbation. 

The  perturbation  expansion  for  the  localized  functions  and 
the  one-band  effective  Hamiltonian  in  this  work  is  of  quite  general 
character  and  can  be  applied  to  a variety  of  problems,  including  the 
impurity  problem,  surfaces,  superlattices  and  external  fields^  ^ . The 
applicability  of  the  approach  is  limited  to  crystals  in  which  band 
touching  does  not  appear  and  there  is  at  least  one  isolated  band. 

Good  examples  of  this  kind  of  crystals  are  the  alkali  halides,  the 
rare  gas  solids  and  the  II-VI  semiconductor  compounds . 

The  idea  of  using  Wannier  functions  in  the  problem  of  a 

g 

Bloch  electron  in  a magnetic  field  was  first  given  by  Luttinger  . 

Later  this  idea  was  further  developed  for  deriving  effective  one-band 

9 

Hamiltonians  in  a paper  by  Roth  . 

The  idea  of  modifying  the  Wannier  functions  in  order  to 

adjust  them  to  the  perturbation  of  a magnetic  field  is  put  forward 

in  the  Refs.  (8)  and  (9).  Following  our  perturbative  approach,  the 

9 

Roth  operator  technique  is  used  for  developing  an  orthonormal  set 
of  Wannier  functions  in  the  presence  of  a magnetic  field.  The  newly 
defined  Wannier  functions  are  shown  to  be  useful  in  deriving  one- 
band  effective  Hamiltonians  for  a Bloch  electron  in  a magnetic  field 
in  an  entirely  elementary  way.  These  functions  are  also  used  to 
prove  that  the  eigenfunction  for  a Bloch  electron  in  a magnetic 
field  can  be  written  as  a product  of  a known  operator  and  the 
envelope  function  which  satisfies  the  effective  Hamiltonian  equation. 
This  structure  of  the  eigenfunction  serves  as  a basis  for  defining 
multiband  operators  in  the  space  of  the  envelope  functions  for 
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different  physical  quantities,  e.g.  the  radius  vector,  the  velocity, 
the  acceleration,  etc.  Different  relations  are  obtained  between 
these  operators  fully  resembling  the  corresponding  relations  in  the 
absence  of  a magnetic  field.  In  particular  an  acceleration  theorem  is 
proven  to  hold  for  the  velocity  of  a Bloch  electron  in  a magnetic 
field.  It  is  also  shown  that  a sum  rule  is  satisfied  by  the  newly 
defined  operators^. 

3.  EQUATION  FOR  LOCALIZED  STATES  IN  PERTURBED  CRYSTALS  (REF.  (6)  ) 

Let  us  start  with  a number  of  remarks  about  the  kq  representa- 
tion and  the  problem  of  localized  states  in  perfect  crystals. 

The  connections  between  a wave  function  i(/(r)  in  configuration 
space  and  the  corresponding  wave  function  C(k,q)  in  the  kq  representa- 
tion are  as  follows’'"'*'  : 


/ v \ 1^2  f 

?>  = to)  / d^c<I•^>• 


u>  -(to 

\ (2tt)3/ 


exp(ik‘R)tKq  - R), 


where  V is  the  volume  of  a unit  cell  in  the  Bravais  lattice,  the 
o 

integration  in  (1)  is  over  the  Brillouin  zone,  and  the  summation  in 
(2)  is  over  all  the  vectors  R of  the  Bravais  lattice.  From  (2)  it 
follows  that  a Bloch  function  Cnj^(lt,q)  in  the  kq  representation  is 

)1'2  *nkWZ.S<J  <« 


i yyjjW'lT1 
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where  , (q)  is  the  Bloch  function  in  the  coordinate  representation 
-*-n  -*■ 

(with  r replaced  by  q),  kg  is  the  Bloch  quasimomentum,  and  the 
summation  is  over  all  the  vectors  K of  the  reciprocal  lattice.  By 
definition,  a Wannier  function  a'nm(k,q)  in  the  kq  representation 
belonging  to  the  nth  band  and  located  on  site  of  the  Bravais  lattice 


/ V f 


= exp(-ik.*Rm)^nk(q)  (4) 

This  is  a very  simple  connection  between  a Wannier  function  in 
the  kq  representation  and  a Bloch  function.  The  orthogonality  of 
the  Wannier  functions  (4)  is  seen  immediately 


dkdqa^k.qja'^,  (k,q) 


= jrexp[ik-(Rm  - Rm,)]^(q)yrn,k(q)d^dq 


= 6 ,6  , 

nn  mm' 


For  arriving  at  (5),  we  used  the  following  relation  that  Bloch 
functions  satisfy  : 


m li 


,k(q)dq  =snn. 


In  what  follows  we  shall  prefer  to  work  with  the  periodic  part 
U(k,q)  of  the  wave  function  C(fc,q)  : 


and  are  the  energy  bands  (kg  being  the  conserved  quasimomentum) . 


The  periodic  part  of  the  Wannier  functions  a^^k.q)  is  therefore 


anm(k’<l) 


exp(-k*Rffi)unk(q) 


(10) 


In  what  follows  anm(k,q)  [and  not  the  a^Ck.q)  in  (4)]  will  be 
called  the  Wannier  functions.  It  is  obvious  that  the  Wannier  functions 
(10)  satisfy  in  the  kq  representation  the  Bloch  equation 


( 


t-lQ/3q)  + k]' 
2m 


V(q)) 


anm(k’<) 


en<k>anm<k^>  <U> 





f 


- 7 - 

It  is  of  interest  to  compare  Eq.  (8)  for  Bloch  functions  with  Eq. 

(11)  for  Wannier  functions.  In  Eq.  (8),  en(kjj)  is  a constant 
independent  of  the  kq  coordinates  while  on  the  right-hand  side  of 
Eq.  (11)  we  have  e^ft),  a function  of  ft.  We  see  that  the  only 
difference  between  the  equation  defining  Bloch  functions  [Eq. 

(8)]  and  the  equation  defining  localized  functions  [Eq.  (11)]  is  in 
the  factor  multiplying  the  function  on  the  right-hand  side.  This 
fact  will  be  used  in  postulating  an  equation  for  localized  states 
in  perturbed  crystals. 

• 4 

The  Schrodinger  equation  for  a Bloch  electron  in  a constant 
magnetic  field  ft  and  perturbation  v(r)  in  the  kq  representation  is^ 

( [-l3/3q  + ft  +(e/2c)ftxi3/3Tc]2  + v(+}  + v/1_ijL(^^)  „ 

V 2m  V W) 

= eu(k,q)  (12) 

Equation  (12)  is  the  eigenvalue  equation  (e  is  a constant)  for  an 
electron  in  a perturbed  crystal.  The  solutions  U(k,q)  are  eigenstates 
of  the  problem.  In  order  to  obtain  an  equation  for  localized  states 
we  shall  use  the  analogy  with  equations  (8)  and  (11).  Equation 
(8)  is  an  eigenstate  equation  [with  a constant  en(kg)  multiplying 
the  function  on  the  right-hand  side],  while  Eq.  (11)  defines  localized 
states  anm(ft,q)  [the  energy  en(k)  multiplies  on  the  right-hand 

side].  Having  this  in  mind  we  shall  postulate  in  correspondence  with 
Eq.  (12)  the  following  equation  for  localized  functions  in  a perturbed 
crystal  : 

( C-i3/3q  + k 4-  (e/2c)ftxi9/9k]2  + 

V.  2m 


+V(1^|))Anm(^^)  = V^^nm^ 


(13) 
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where  A (k,q)  are  the  perturbed  localized  functions  and 
nm 


En(^,Rm)  are  the  effective  band  energies.  The  explicit  dependence 
on  the  site  R^  appears  in  the  equation  if  one  follows  the  structure 
of  localized  functions  (10)  in  the  kq-representation 


A«(k.q>  - exp(-ik-Rm)Unm(k,q) 


with  U^Ck.q)  satisfying  the  equation 


(14) 


( 


[-i3/3q  + k + (e/2c)Hx(R  + i3/3k)]2 

PI 

2m 


+ V (q)  + 


v(i—  + R U 

k di  nm 


(k,q) 


En(k’Rm>°nm<k’q> 


(15) 


In  the  absence  of  perturbations  (H  = 0,  v = 0)  Eq.  (15)  goes  over  into 

Eq.  (11)  for  localized  states  of  a perfect  crystal.  In  the  latter 

case  the  localized  functions  U (k,q)  and  the  effective  band  energies 
+ + -►  nm 
E (k,R  ) are  site  R independent 
n m m 


Unm(k,q)  = unk(q) 


(16) 


E (k,R  ) = e (k) 
nv  ’ nv  ' 


(17) 


This  is  no  longer  so  when  a perturbation  is  present.  As  is  seen 
from  Eq.  (15),  both  the  localized  functions  UnTn(k,q)  and  the  effective 
band  energies  E (k,R  ) will,  in  general,  be  R dependent  (an  exception 
is  the  homogeneous  electric  field  case).  One  should  expect  Eq.  (15) 
to  lead  to  localized  states  because  the  solution  of  the  unperturbed 
equation  [Eq.  (11)]  are  localized  and  they  should  not  be  strongly  affected 
by  the  perturbations.  It  is  therefore  possible,  in  principle,  to  solve 


I 
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Eq.  (15)  by  a perturbation  procedure  starting  with  the  unperturbed 

equation  (11) . This  Is  the  main  advantage  of  the  localized 

equation  (15)  over  the  eigen-value  equation  (12).  The  eigen- 

functions  U(k,q)  of  the  latter  are  very  different  from  the  eigen- 

functions  U , (k,q)  of  the  unperturbed  problem  (8)  and  a perturba- 
T“Cb 

tion  procedure  can,  in  general,  not  be  applied  to  Eq.  (12). 

For  comparison  reasons  let  us  write  Eqs.  (11)  and  (13)  for 
localized  states  in  the  regular  r representation.  The  left-hand  sides 
of  Eqs.  (11)  and  (13)  will  just  be  the  corresponding  Hamiltonians 

-f 

applied  to  the  function  in  the  r representation.  In  order  to  obtain 

the  right-hand  side  we  expand  the  band  energies  e (k)  and  E (k,R  ) in 

n n m 

Fourier  series  and  use  formula  (1) . Equations  (11)  and  (13)  in  the  r 
representation  will  be 


( 4 + v(?) 


>) 


Vr 


->■ 

- R ) 
m ' 


+ S8>' 


( 


[p  + (e/2c)Hxr]2 
2m 


+ V(r)  + v(r) 


) 


A (r  - R ) 
nmv  m 


y H (R  R )A  (r  - R + RJ, 
nx  s,  m nm  m s ’ 


(19) 


where 


W 3 

n s (27 r)J 


( -iic*R 
)e  8£n 


(k)dk 


H (R  R ) 
n s,  m 


(2tt)  ' 


j e ik  RsEn(k,Rm)dk 


(20) 


(21) 


and  a(r-R),A  (r-Rn)  are  the  localized  functions  for  the 
n m nm  jg 

perfect  and  perturbed  crystal,  correspondingly.  Equation  (18)  is  well 
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known  and  is  the  equation  for  Wannier  functions  aQ(r)  In  a perfect 
crystal.  In  this  case  the  shape  of  the  function  afl(r  - Rn)  does 
not  vary  with  the  site  and  all  the  Wannier  functions  for  a given 
band  are  derived  from  a single  function  a (r).  In  the  perturbed 
crystal  the  shape  of  the  localized  function  Ar||n  (r  - R^)  depends  also 
on  the  site 

/ V x1/2  f 

W?  ' ■ (7-T3)  J «*P(-ik'Vc0,(k.')  <22> 

where  C (k,r)  * exp(ik«r)U  (k,r)  [see  Eq.  (1)].  In  the  regular 
ran  mu 

^representation,  Eqs.  (18)  and  (19)  appear  to  be  more  complicated 
than  Eqs.  (11)  and  (13)  and  we  find  it  more  convenient  to  work  In 
the  kq  representation. 

4.  THE  IMPURITY  PROBLEM 

In  the  presence  of  an  impurity  perturbation  v (H  ■ 0),  Eq. 
(15)  for  localized  states  becomes 


v(?)  + v (i-i  + 


' ‘.W.'VM  (23) 

We  shall  solve  Eq.  (23)  by  a perturbation  procedure  assuming  that  the 
unperturbed  bands  £ (k)  are  well  separated.  This  will  be  needed  for 
defining  a parameter  of  smallness  for  the  perturbation  theory.  Any 
function  U(k,q)  can  be  expanded  in  the  complete  set  of  functions  (10) , 


U(£,q) 


% «*P(-lk-R.)unk(5) 


• VE)u»k* 
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Hr 

r ■ 

r*  ■***—*«- 


u 


which  is  a well  known  expansion  in  the  kq  representation  . For 


the  solutions  U (k,q)  of  Eq.  (23),  expansion  (24)  will  be 

luu 


%<*.«>  ’ 21 

l 


(25) 


In  this  case  the  expansion  coefficients  are  site  Rm  dependent.  We 
substitute  expansion  (25)  into  Eq.  (23),  multiply  it  by  u*fc(q)  from 
the  left,  and  integrate  over  q.  We  arrive  at  the  following  equation 
for  the  coefficients  Bnl(k*V  ! 


(26) 


where 


^ us*(q)v  (i-|  + Rm)utk(q)d^ 


(27) 


Expression  (27)  defines  an  operator  which  is  to  be  applied  to  the 
coefficients  Bn^(^»\)»  80,1  the  integrand  from  point  of  view  of  the 
quasimomentum  t is  a product  of  three  operators.  Equation  (26)  can 
be  solved  by  perturbation  theory.  Let  us  introduce  the  following 
notations  : 


BM<k> 


+ Bns>(k’*.>  + + ••• 

(28) 


En(t,^)  - en(k)  + (£,$b)  + 42)(£.$b)  + *•*  <29> 


U 
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where  the  superscript  denotes  the  order  of  perturbation.  Up  to 
second  order  in  the  perturbation  we  have  (we  apply  perturbation 
theory  to  band  n) 

B^°*(k,R  ) = 6 
ns  * nr  ns 


Bn^  (k,  R)  = 8g(  -m)..  , s * n.  B^k.R^)  - 0 (31) 

ns  m e (k)  - e (k)  1111  m 

11  S 


B(2)(k,R  ) 
ns  ' * nr 


vst(k,i3/35  + 


*n00  - e.<k> 


vnn(R-R.»)v6..(R’V 

- es<f>2 


8 + n 


_ 1 Vtn(t-R„> 

2 2-  [E  (I,  . e (!)]2 
Wn  " ‘ 


V <k'V  ■ 'W<k-Rm)- 


^2)  <*•*»>  - z + K)  *2’  <*•*■>  <35> 

3k 

In  the  above  formulas  in  addition  to  the  matrix  operator  (27)  also  the 
matrix  v£8(k,Rm)  appears.  The  latter  is  a special  case  of  the 
operator  (27)  when  it  no  longer  operates  on  a function  of  £ [the 


L 


operator  v(i3/3k  + Rm)  in  the  integrand  of  (27)  operates  in  this 
case  directly  on  the  function 

vs*<k»Rm>  “ /[_  exp(-ik*R)  x dr a* ( r-R) v ( r+Rm) a^  (r  ) 

R 

- 21  exp(ik*R)v8Jl(Rm+R,Rm)  (36) 

R 

with 

v^R.R’)  = Ja^-RM^a^-R’M?  (37) 

For  results  (30)- (35)  to  have  the  meaning  of  perturbation  corrections 

we  have  to  assume  that 

n = vsn(^V/ten(^)-£s(^)]  <<  1*  (38) 

where  0 is  in  some  sense  the  expansion  parameter  of  the  perturbation 
theory. 

When  (38)  holds  the  localized  functions  An|n(k,q)  for  a crystal 
with  an  impurity  up  to  second  order  in  perturbation  theory  will  be 

*„<£.?>  - e-'PC-lS'V  *(fl  + B“>(S.*,))»nk(?) 

+ Z + 

s^n 

where  the  coefficients  B^^.R^)  are  given  by  (31)-(33). 

The  functions  (39)  by  the  construction  of  perturbation 
theory  are  normalized  to  second  order  in  the  perturbation 
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.<{.?)  1 2dM?=l  + B<2)(I.Rm,  ♦ »«>*(£,  ^ 


+ Z-  - 1 0 

s In 

The  last  equality  in  (40)  Is  a consequence  of  (31)  and  (33). 

By  using  these  results  a one  hand  effective  Hamiltonian 

can  be  derived  for  the  impurity  problem.  The  derivation  will  be 

based  on  the  localized  functions  A (Icq)  [Eq.  (39)]  and  the 

effective  one  band  energies  (39)  and  (35).  Following  some  simple 

physical  assumptions  of  Ref.  (6),  the  localized  functions  A (tq) 

nm 

assume  the  following  form  to  first  order  of  the  perturbation 


* exp(-ik-Rm)  x^unk(q)  + 
y vns  -*  ) 

Z-  -T—  V0»J 

j i ns  " 


The  structure  of  the  localized  functions  (41)  is  very  simple  and  they 
can  easily  be  written  in  the  r representation 

Anm(*  ' *m>  " an(*  ‘ V + 21  a8<*  ’ *„>  <*2> 

sin  tns 

Here,  Anm(r  - R^)  are  the  localized  functions  for  a crystal  with  an 

Impurity.  It  is  seen  that  A (r  - R ) depends  not  only  on  the 

.+  ““  m ^ + 

difference  r - as  in  the  unperturbed  Wannler  functions  aj  - V 
but  also  explicitly  on  the  site  R^.  As  should  be  expected,  the 
stronger  the  perturbation  is  on  site  R , the  more  will  the  perturbed 
function  A l L(r  - Rffl)  differ  from  the  unperturbed  one  ®n(r  - on  this 
particular  site  R . 

ID 

The  assumptions  simplify  also  considerably  expressions  (34) 
and  (35)  for  the  effective  energies.  To  second  order  in  perturbation 
theory  we  have 
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V • 'JK?  ■ 


t>  - Z 


Iv  (R  ) I : 
ns  m 


s^n  na 


Let  us  now  construct  a one-band  effective  Hamiltonian  for  the 


impurity  problem.  For  this  purpose  we  expand  the  solution  U(k,  q)  of 


Eq.  (12)  (with  H = 0)  in  the  localized  function  A^(k,  q) , 


»»•<>-£  YVVk>  ’>• 


where  the  expansion  coefficients  FIR  ) depend  on  the  band  index  Jt  and 
-*■  IP 


the  site  R . Substituting  (45)  into  (12)  we  have 
m 


2L  v1*  VYVVk-  <>  ‘ F*(WE-  ’>  • 


where 


-+  -►  "►  lv0e(R  ) I 

h(K  V ■ Yk)  + v.  <V  + Z. 

F * p sH  is 


Up  to  the  order  of  U the  one-band  effective  Hamiltonian  equation  for  the 
impurity  problem  will  be 


h (R  - R )F  (R  ) 
n p m n p 


+ (v  (R  ) + 
nn  m 


Iv  (R  ) I2  ->  + 

-T  )F.(»J  - eFn(R_ ) . 


syn  ns 


Equation  (49)  has  a correction  term 


Iv  (R  ) I 
ns  m 


s^n  ns 


to  the  well-known  one-band  classical  equations  for  localized  impurities  of 

» j 1 l r * . « j t ^ . _ i e r j 


Roster  and  Slater  and  for  shallow  impurities  in  the  effective-mass 
approximation.  Equation  (48)  reproduces  both  mentioned  equations  when 


J 
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the  correction  term  is  neglected.  However,  the  second-order  perturbation 
term  in  Eq.  (48)  is  of  very  simple  structure  and  can  lead  to  significant 
contributions  in  the  one  band  impurity  potential.  The  significance  of 
this  term  follows  from  the  fact  that  while  being  a one-band  band  term 
it  is  influenced  by  the  explicit  band  structure  of  the  solid.  This  can  in 
particular  be  easily  seen  in  solving  the  Koster-Slater  localized  impurity 
problem.  Thus  by  assuming  that  only  two  bands  Interact  in  Eq.  (48) 

(band  s influences  band  n)  and  that  the  impurity  is  localized  at  the 

origin,  we  find  that  the  effective  impurity  potential  at  R >0  is 

n 


V 

nn 


(0)  + 


|v  (°)|2/ena 
'ns  1 ns 


(49) 


This  potential  depends  explicitly  on  the  band  gap  Cqs  and  the  interband 

matrix  elements  v (0) . The  second-order  correction  term  in  (49)  can  be 
ns 

large  and  is  restricted  only  by  the  applicability  of  perturbation  theory. 
Both  the  binding  energies  of  the  impurity  and  scattering  cross  sections 
will  be  modified  by  the  band  structure. 

Similar  remarks  can,  in  principle,  be  made  about  the  significance 
of  the  correction  term  in  Eq.  (48)  for  the  shallow  impurity.  In  this  case, 
however,  one  should, in  general,  expect  that  for  gentle  potentials  v(r) 
the  second  order  term  in  (48)  will  be  very  small.  One  could,  however, 
imagine  cases  when  vs£^)  is  not  negligibly  small.  This  is  only  a guess 
and  it  should  be  checked  on  real  problems. 


5.  BLOCH  ELECTRONS  IN  A MAGNETIC  FIELD  (REF  . (10)) 

In  the  presence  of  a homogeneous  magnetic  field  H the  equation  for 
localized  states  becomes 

(30) 


<-*  4 + fe Hxi 


v 


2m 


»<*>!  *„<*•»»>  A„&> 
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I 


L 


where  A (kq)  Is  the  perturbed  localized  function  for  band  n and  site  R 
nm^  m 

and  E (1c, R ) Is  the  site  dependent  effective  energy.  In  the  absence  of 
n m 


the  magnetic  f ield.Eq. (50)  goes  over  into  Eq.  (11)  and  correspondingly 


A A 


A (kq)  go  into  a_  (kq)  (Rel.  (10))  and  the  E (k,R  ) become  the 
nm  nm  ^ n m 


unperturbed  energies  c (k) . One  can  look  for  A In  the  form 


nm 


Anm(kq)  = exP(-ik*V  Unm(kq) 


m nm 


(51) 


where  Unni(kq)  satisfies  the  equation' 


(-i  + k + 

[ 19 


nm 


En(k,R  ) U (kq) 
n m nm 


(52) 


The  Hamiltonian  in  (52)  depends  on  the  combination 


k = k + tt-  HxR 
m 2c  m 


(53) 


which  means  that  (51)  can  be  rewritten  as 


Anm(kq)  = exP<-ik'V  Unk  <9) 

m 

* » 

and  E (k,R  ) in  (52)  will  be 
n m 


(54) 


E (k,R  ) 
n m 


E_(k  ) 


n m 


(55) 


The  results  (50)  - (54)  can  be  written  in  a more  compact  form  by 

A ^ A g Jfe  0 

using  the  notation  A(X)  for  a function  of  the  operator  x “ k + Yc  3? * 

A(x)  is  defined  as  follows 


A(X)  - / dA  F(A)  exp{i(k  + ^ H*i  |gJ-A} 


(56) 
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where  F(X)  is  the  Fourier  transform  of  A(k) . In  Ref.  (9)  a simple  multi' 
plication  rule  was  shown  to  hold  for  a product  of  two  operators 


B(X) 


where 


C(k)  - A(k)  B(k)  = ih  0 - 


a8  3ka  3ke 


h 92A(k)  32B(k)  L 

' 2 naf5  a’8'  ak^,  Sk^, 


a8  a8y  2c 

In  (59)  e * is  the  antisymmetric  unit  tensor  and  in  (58)  and  (59)  there 
apy 

arm  summations  over  repeated  indices.  Rel.  (57)  together  with  (58)  and 
(59)  give  what  is  called  the  commutator  expansion  in  powers  of  the  magnetic 
for  a product  of  two  functions  of  x* 


By  using  notation  (58)  the  localized  functions  (54)  can  be  written 


as  follows: 


A (kq)  = U (Xq)  exp(-ik*R  ) (60) 

nm  n m 

»»  * * 

where  un(Xq)  operates  on  exp(-ik*Rm).  Similarly,  Eq.  (50)  can  be  rewritten  as 
H(xq)  Un(xq)  expC-ik-RJ  = (En(k)  U/^k-X  exp(-ik*Ra)  (61) 


where  H(kq)  is  the  Bloch  Hamiltonian.  On  both  sides  of  Eq.  (61)  there  are 

I * 

operators  which  are  applied  to  the  function  exp(-ik*R  ).  Without  the  latter 

n 

Eq.  (61)  is  an  operator  equation.  For  reasons  that  will  become  clear  below  It 
is  convenient  to  write  this  operator  equation  in  a somewhat  different  form 


H<Xq)  UQ<Xq)  " VXq)  Hn(X) 
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On  Che  right  hand  side  of  Eq.  (62)  the  operators  are  rearranged  in  comparison 


with  Eq.  (61)  and  this  is  the  reason  for  writing  H (k)  Instead  of  E (k). 

n n 


Eq.  (62)  will  be  the  main  equation  for  defining Cannier -inactions  and  one- 
band  effective  Hamiltonians -far  •*  Bloeh  electron  in  a -magnetic  field. 


One  can  now  prove  that  the  functions  A ■ (kq)  in  (60)  are  the 

tun 


Wannler  functions  of  the  problem. 


Before  turning  to  the  solution  of  Eq.  (62)  let  us  show  its  con- 
nection with  the  eigenvalue  problem  for  a Bloch  electron  in  a magnetic  field. 
The  SchrBdinger  equation  for  the  latter  in  the  kq-representation  can  be 


written  as  follows 


H(Xq)  U(kq) 


E U(kq) 


where  H(kq)  is  the  Hamiltonian  and  H(Xq)  is  obtained  from  H(kq)  by 
replacing  in  it  k by  k + -5—  Hxi  . It  can  be  assumed  that  the  Warn 


Wannler 


functions  A (kq)  in  (61)  form  a complete  set  and  the  following  expansion 

nm 


holds  for  the  eigenfunction  U(kq) 


U(kq)  * X F (R  ) A (kq)  - 2-  U (Xq)  Bn(k) 
n m nm  n n 

nm  n 


where  B (k)  and  F (R  ) are  Fourier  transforms  of  one  another.  Let  us 
n n m 


Substitute  Exp.  (64)  and  (63)  and  use  Eq.  (62).  We  find 


Z Un(Xq)  H (X)  Bn(k)  = eX.Un(Xq)  B^k) 
n n 


By  multiplying  Eq.  (65)  by  Uff(xq)  from  the  left,  and  by  integrating  it  over  q 


one  finds 


Ha(X)  Bjl(k)  - e B^k) 


which  is  a one-band  effective  Hamiltonian  equation  for  a Bloch  electron  in  a 
magnetic  field.  Eq.  (66)  shows  that  H (k)  in  the  equation  for  localised 
states  (Eq.  (62))  is  the  effective  Hamiltonian  of  the  problem.  This  means  that 
by  knowing  the  solutions  of  Eq.  (62)  one  knows  directly  the  one-band 


r- 
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effective  Hamiltonian.  In  addition,  from  Eq.  (66)  together  with  Exp.  (64) 

one  can  find  the  eigenfunctions  for  a Bloch  electron  in  a magnetic  field. 

Since  the  operators  U (Xq)  diagonalize  the  Hamiltonian,  only  one  term 

n 

should  be  kept  in  the  Exp.  (64)  and  the  eigenfunction  U(kq)  corresponding 
to  the  band  n becomes 


U(k?) 


Un(Xq}  Bn(k) 


(67) 


In  the  kq-representation  the  eigenfunction  U(kq)  exhibits  therefore  a 
separation  of  the  envelope  function  B (k)  and  the  crystalline  operator 


n 


Un(xq).  Rel.  (67)  shows  that  in  the  kq-representation  the  separation  of 


the  eigenfunction  into  a product  of  a known  crystalline  operator  and  an 
envelope  function  is  a general  one  band  feature  without  any  additional 
assumptions.  This  general  feature  of  the  kq-representation  is  a consequence 
of  the  fact  that  k and  q are  the  symmetric  coordinates  in  an  ideal  crystal. 


The  solutions  of  the  main  equation  (Eq.  (62))  can  be  found  as  a 
power  series  of  the  magnetic  field  in  the  framework  of  the  commutator  expansion 
(57)  - (59).  It  was  already  shown  that  these  solutions  give  the  one  band 


effective  Hamiltonian  and  the  Wannier  functions  in  the  presence  of  a 

t 

magnetic  field.  For  solving  Eq.  (62)  multiply  it  by  IL  (Xq)  from  the  left 
and  integrate  over  q 


(2tt) 


3 V(X)  = 1 VXq)  H(Xq)  Un(Xq)dq  = 6Jtn 


(2tt) 


3 VX> 


(68) 


The  next  step  is  to  expand  the  function  U (kq)  in  Wannier  functions  of 


the  ideal  crystal 


un(kq)  - 2-  Bns(^)  u«k(<l> 


sk' 


(69) 


The  unknowns  are  now  the  effective  Hamiltonian  H (k)  and  the  coefficients 
A n 


B (k).  One  can  look  for  them  in  the  form  of  a perturbation  series 

OB 


H (k)  - e(o)(k)  + e(1) (k)  + e(2)(k)  + ... 
n n n n 


^ 111  (2)  ^ 

B (k)  - 6 + BV  ; (k)  + BV  ' (k)  + ... 

ns  ns  ns  ns 


(70) 

(71) 


r 
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where  the  superscripts  denote  different  orders  in  the  magnetic  field 


« 

t 

according  to  the  commutator  expansion  (57)  - (59).  It  can  be  shown  (Ref.  (10)) 
that  the  effective  Hamiltonian  (Rel.  (70))  and  the  coefficients  (Rel.  (71)) 
can  be  determined  to  any  order  of  the  magnetic  field  from  Eq.  (68)  together 
with  the  orthogonality  relation.  By  solving  these  equations  it  is  shown  with 
what  ease  one  can  reproduce  the  well  known  results  of  constructing  one-band 
Hamiltonians  for  a Bloch  electron  without  the  need  of  using  diagonalization 
procedures,  Thus,  the  zero  order  effective  Hamiltonian  is 

40)<*>  ■ 

(72) 

while  the  first  order  results  are 

e<1}  (k)  = V^}  (k) 

(73) 

. 

00  = - J (k) 

(74) 

1 % 

while  for  9 + n one  finds 

.(us,  _ C<*> 

nA  en(k)  - 

(75) 

The  effective  Hamiltonian  (73)  is  identical  with  the  one  given  by  Formula  (73) 
of  Ref.  (9).  It  is  interesting  that  both  the  effective  Hamiltonian  (Ref.  (73)) 
and  the  coefficient  (Rel.  (75))  to  the  first  order  in  the  magnetic  field  are 
determined  by  a unique  non-Hermitian  matrix.  The  expressions  (73)  and  (75) 
look  like  first  order  perturbation  formulas. 

• 

With  the  same  ease  one  finds 

(see  Ref.  (10)): 

the  second  order  terms  in  (70)  and  (71) 

• 

• 

e^2)(k)  - G<2)(k)  - e (k)  N(2)(k) 

n nn  n nn 

1 (k)  1 ^ 

+ T_  _ „<!>£)  V<D£) 

tr„  *„«>  - s*>  *- <k)  v-» <k) 

(76) 

A 
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’ ‘ 1 <">  + l»"«l2  + ^ at  Z.  »£>»>  xa„<‘)  * ».c.) 

S p 8 

(77) 


By  using  the  latter  expressions  one  can  show  that  the  effective  Hamiltonian 
0 in  (76)  is  identical  with 

Similarly,  for  % + n we  have 


(2) 

en  '^)  In  (76)  is  identical  with  the  one  given  by  Formula  (74)  of  Ref.  (9). 


- en<k>  <k>  + 

+ X-B.(i>*(k)  B^OO  (e0(k)  - e. (k) ) + (B^ab  Vfl(J}(k)  + h.c.)  + 


+ Ho$  9k~  (i  Bn£  (k)  v6U(k)  + h,c,)  + ha6  3k“  X Bns)(k)  *8*8 (k)  (78) 

01  Ct  8 

It  is  therefore  seen  that  with  very  little  algebra  one  finds  from 
the  Eq.  (44)  both  the  effective  one-band  Hamiltonian  and  the  coefficients 
Bni  (k)  . The  latter  determine  the  Wannier  functions  in  a magnetic  field 
(Reis.  (61)  and  (69)).  Thus,  to  the  zero  order  of  the  magnetic  field  one  has 

f*a°i  - V 


Anm(k<») 


unk  («>  e*P<-ifc*Rm> 


(79) 


m 


where  k is  given  in  (53).  The  zero  order  Wannier  functions  in  the  presence 

of  a magnetic  field  (Rel.  (79))  are  very  similar  to  the  corrmaponding  functions 

in  an  ideal  crystal.  The  former  are  obtained  from  the  latter  by  replacing  k 

in  the  periodic  function  u . (q)  by  k . The  functions  (79)  are  modified 

nx  m 

Wannier  functions  for  the  magnetic  field  problem  and  they  were  already  used  for 

8 q 

describing  the  motion  of  a Bloch  electron  in  a magnetic  field.  ' Similarly, 
one  obtains  the  Wannier  functions  to  higher  order,  of  the  magnetic  field. 
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In  summary,  the  use  of  localized  functions  that  correspond  to  non- 
ideal crystals  enables  one  to  take  directly  into  account  multiband  effects 
in  the  impurity  problem  and  in  the  problem  of  external  fields.  Wannier 
' functions  in  ideal  crystals  have  been  known  as  an  extremely  useful  tool  in 

> general  theoretical  proofs  and  model  calculations  in  solid  state  physics. 

' > The  work  reviewed  in  this  report  shows  that  a similar  role  is  played  by 

perturbation  adapted  localized  functions  in  imperfect  crystals. 
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